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Abstract 

We study the decay f2 _ — ► H~7r + 7r~ in heavy-baryon chiral perturbation theory. At leading 
order, the decay is completely dominated by the H*°(1530) intermediate state, and the predicted 
rate and H~7r + -mass distribution are in conflict with currently available data. It is possible to 
resolve this conflict by considering additional contributions at next-to-leading order. 
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I. INTRODUCTION 



It was suggested many years ago that the decay Sl~ — > 2~7r + 7r~ should be dominated by 
the S*°(1530) intermediate state [l|, 0]. Under this assumption, the current Particle Data 
Group [3| branching ratio for Q~ — > H*°7r~ has been deduced from the measurement of 
£>(fi~ — > 2~7r + 7r~) [4]. More recently, the HyperCP collaboration has reported a preliminary 
measurement of Q~ — > 2~7r + 7r~ that is very surprising in that the distribution of the 5~7r + 
invariant-mass apparently shows no evidence for the 5*° (1530) dominance 0. 

Motivated by this result, we revisit the calculation of the rate for this decay mode using 
heavy-baryon chiral perturbation theory (HB^PT). We first present a leading-order cal- 
culation that reproduces the expectation that the decay is completely dominated by the 
S*°(1530) intermediate state. 

We next explore whether higher-order contributions can reconcile the calculation with 
the preliminary HyperCP result. To this end, we consider the effect of next-to-leading-order 
diagrams, which occur at tree level. 



II. LEADING-ORDER CALCULATION 



The amplitude for Q 
approach as 



n TT TT 



-u 



H (p s ) 7r + (j» + ) tv (p_) can be written in the heavy-baryon 



a (A + pt + A_p»+ 2B + S v ■p_p» + + 2B_ S v -p +P t) u n ^ , (1) 



where A± and B± are independent form-factors and S v is the spin operator. The most 
general form of the amplitude has eight independent form-factors |l( > an d we have included 
here only the ones that receive contributions from the leading-order and next-to-leading- 
order diagrams that we consider. The partial decay width resulting from the amplitude 
above is 



where mSL + 



1 



C TT TT 



32 (2vrm n )' 



\M(Q- 



C TT TT 



'• dml 



dm E-7T- , 



(2) 



(Ps + P±) 2 an ^ 



\M{n- 



^~TT^TT~ 



{ \A + \ 2 p 2 + + \A_\ 2 p 2 _ + 2Re(A* + A_) p + p 



\A 
+ 



\B, 



f \B_\ 2 + Re(B* + B_) 
+ Re(BlB_) (p + -p_) 2 } , 



„2 2 
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with p ± denoting the three-momenta of the pions in the Q rest frame. 

The chiral Lagrangian describing the interactions of the lowest-lying mesons and baryons 
is written down in terms of the lightest meson-octet, baryon-octet, and baryon-decuplet 
fields The meson and baryon octets are collected into 3x3 matrices cp and B, 
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respectively, and the decuplet fields are represented by the Rarita-Schwinger tensor T£ bc , 
which is completely symmetric in its SU(3) indices (a, b, c). The octet mesons enter through 
the exponential S = £ 2 = exp(iip/ f), where / = f n = 92.4 MeV is the pion-decay constant. 

In the heavy-baryon formalism [sj], the baryons in the chiral Lagrangian are described by 
velocity-dependent fields, B v and T£. For the strong interactions, the Lagrangian at lowest 
order in the derivative and m s expansions is given by 

C, = (4i^(c>^+[V M ,£j))+2^ 

- T£ iv ■ VT V)l + Am T$T vll + C (T V ^B V + B V A^) + 2H T»S V ■ AT Vfl 

+ b D (B v {M + , B v })+b F (B v [M + , B v ] > + cT»M + T vtl (4) 

where only the relevant terms are shown, (■ • ■ ) = Tr(- • • ) in flavor-SU(3) space, Am de- 
notes the mass difference between the decuplet and octet baryons in the chiral limit, V M = 

and M+ = ftjvf£t + £M f £, with M = diag(m, m,mj = diag(m 2 , ml, 2m 2 K -rnQ/(2B ) 
in the isospin-symmetric limit m u = m d = m. The constants D, F, C, Tl, Bo, bp^, and c 
are free parameters which can be extracted from data. 

As is well known, the weak interactions responsible for hyperon nonleptonic decays 
are described by a (AS*! = 1 Hamiltonian that transforms as (8l, 1r) © (27l, 1r) under 
SU(3)lxSU(3)r rotations. It is also known empirically that the octet term dominates the 
27-plet term. We therefore assume in what follows that the decays are completely charac- 
terized by the (8 l ,1r), |A/| = 1/2 interactions. The leading-order chiral Lagrangian for 



such interactions is 

C w = h D (B v {?hZ,B v }) + h F (B v [?hZ,B v ]) + hof£tfhZT vll + H.c. , (5) 

where h is a 3x3 matrix having elements hki = 8^31 and the parameters ho,F,c can be 
fixed from two-body hyperon nonleptonic decays. 

From C w together with C s , we can derive the O(p ) diagrams displayed in Fig. [TJ They 
provide the leading-order contributions to the A± and B± form factors in Eq. ([TJ, namely 

4 0) = 9 , +Ch ° r , (6a) 

A { ° } = , (6b) 

B (o) = -CHh c 

+ 18 / 2 (m n -m E ,)(E E + E + - m H ») ' 1 ' 

B (o) _ ~C(D- F) h c 



C7ih c 

+ 27 f 2 (m n - mj) (E 3 + E + - m s .) ' ^ ' 
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FIG. 1: Diagrams contributing to fi 
sents a strong vertex from £ s in Eq. i 



-> S 7r + 7r at leading order in xPT. Each solid dot repre- 
) , and each square a weak vertex from £ w in Eq. ([5]) . 



where m- 



Numerically, to evaluate the decay rates resulting from the form factors above, we employ 
the tree-level values of the strong and weak parameters. Specifically, 



D 



0.80 



0.46 



\C\ 



1.7 



(7) 



from hyperon semileptonic decays and the strong decays T — > B(p, but a tree-level value 
of Ji is not yet available from data. Since nonrelativistic quark models give 3F = 2D, 
C = —2D, and Ti, = —3D, which are well satisfied by D, F, and C, we adopt 



n 



-2.4 



For the weak parameters, we have 



h c = 3.42 x 10" 



GeV 



(8) 



(9) 



hj^ = — 1.45x 10 -8 GeV, and hp = 3.50x 10 -8 GeV, extracted from a simultaneous tree-level 
fit to the 5*- wave octet-hyperon and P-wave Q~ nonleptonic two-body decays, as hp^ con- 
tribute not only to the octet-hyperon decays, but also to Q~ — > AK, whereas h c contributes 
to Q~ — > AK, S7r [9]. As seen above, he is the only weak parameter in the lowest-order 
contributions to Q~ — 



The resulting branching ratio, 



C IT TT 



5.4 x 10" 



(10) 



is roughly an order of magnitude larger than the preliminary number reported by HyperCP, 
B{Qr -> S-7T+7T-) = [3.6 ± 0.3(stat)] x 10" 4 and also the current PDG value, B(Qr -> 
S _ 7r + 7r _ ) = (4.3lf 4 ) x 10~ 4 In Fig. E^a), we display the corresponding S~7r + invariant- 
mass distribution. As expected, these results are dominated by the H* resonance. Notice 
that the leading-order rate is proportional to \Chc\ 2 so that there is a large parametric 
uncertainty in this prediction. For example, if both C and he were 30% smaller than the 
values we used, the predicted rate would be four times smaller. The general dependence of 
the leading-order branching ratio on \Chc\ is shown in Fig. |2fb). 
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FIG. 2: (a) Distribution of E~tt + invariant-mass in £l~ — > E _ 7r + 7r~ at leading order with parame- 
ter values in Eqs. ©-([8]), and (b) its branching ratio as function of \Chc\ with D—F and Ti values 
in Eqs. © and ©. 



The HyperCP data is not available in a format suitable for direct comparison with our 
result due to detector effects. However, their results indicate that a uniform phase-space 
distribution is a much better fit to the data than a 2*-dominated one In Fig. [3] we 
plot the 777.=-^+ distributions resulting from our leading-order amplitude (solid curve) and 
from assuming a uniform-phase-space decay distribution (dashed curve), both normalized 
to reproduce the central value of HyperCP's result. The structure of the leading-order 
amplitude, from Eq. ([6]), with all the terms being proportional to Che-, is such that the 5* 
resonance is always the dominant feature of the spectrum. This leads us to investigate in the 
next section whether any of the next-to-leading-order corrections can modify the predicted 
spectrum in the direction indicated by experiment. 
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FIG. 3: Distributions of H~7r + invariant-mass in 0,~ — > E~7r + 7r~ obtained from our leading-order 
amplitude (solid curve) and from the assumption of uniform-phase-space decay distribution (dashed 
curve), both normalized to yield B(Q~ — ► H~7r + 7r~) = 3.6 x 10~ 4 . 



III. CALCULATION TO NEXT-TO-LEADING ORDER 



At next-to-leading order, 0(p), there are two types of contributions. The first type of 
contributions is that in which the weak transition occurs only between mesons. To compute 
these contributions, we need the leading-order, 0(p 2 ), strong and weak Lagrangians for 



mesons, which are given respectively by p, [10 



4 = \f{d^d^) + \B Q f{M + ) 



(11a) 



4, = ij 2 {hd^d^) + H.c. 

where the parameter 7 8 is found from K — > tttt data to be 

7 8 = -7.8 x 1CT 8 , 



(lib) 



(12) 



the sign following from various predictions 11 . 

The contributions of the 7 8 term are interesting because the | AS\ = 1 weak transitions in 
the meson sector are larger than naive expectations. In particular, 7 8 is several times larger 
than its naturally expected value (~1 x 10~ 8 ) and therefore could make its contributions 
numerically comparable to the lower-order ones. 

With weak vertices from the 7 8 term alone, plus strong vertices from C s and C s , we 
derive the next-to-leading-order (NLO) diagrams displayed in Fig. HI They provide the 
NLO contributions to the A± and B± form factors in Eq. (flj), namely 
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(i) = ~ C % m l ~ s +- 
f 2 m\ - s. 
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FIG. 4: Diagrams contributing to 17 _ — > H~7r + 7r _ at next-to-leading order in xPT. Each solid dot 
represents a strong vertex from £ s in Eq. @ or C' s in Eq. (jllaj) . and each square a weak vertex 
from C' w in Eq. ()1 lb|) . 

A w = , (13b) 

(13 -CW Tgm^ , . 

+ 3/2 ( m * K - r n*){E s + E + -m B .)' [ ' 



B 



(i) 



/ 2 (m^ - m 2) (E s + E + - m B ) 

+ ?™ 1*< (13d) 

^ 9f 2 {m* c -ml)(E 3 + E + -m B .) ' 1 J 

There is another type of NLO contribution to the amplitudes. It is given by diagrams 
similar to those in Fig. [1] in which one of the vertices is from a NLO Lagrangian. Many of 
the parameters in NLO Lagrangians are not known, and so it is not possible at present to 
include their contributions in a detailed way. For example, the weak Lagrangian at 0(p) 
that generates Q~E*tt and 17~H7r vertices is, as discussed in Appendix |A~1 

+ ^d^ + E°n- + ••• , (14) 

where only the relevant terms are displayed and hns*n, hnE*n, and hnEn contain unknown 
parameters. The vertices occur in diagrams similar to the first one in Fig. [T]with intermediate 
S* and S, yielding the NLO contributions 

i« = r - Ch ^ E ~ , (15a) 

i W = , (15b) 
= 7= r Ch ^ r , (15c) 
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V2f 2 (E 3 + E + -m s ) 3V6f 2 (E E + E + -m H .) 
Numerically, we adopt the parametric variations 

< \h us *„\, \h us *„\, \h nS7T \ < 2xl0" 8 , (16) 

where the upper limit is the expectation from naive dimensional analysis. 

As mentioned above, there are additional NLO contributions that are not included in our 
calculation because they depend on more unknown parameters. We can still estimate the 
uncertainty in our results arising from those terms by allowing the LO parameters to vary 
between their value as obtained from tree-level fits and their value as obtained from one-loop 
fits. For our numerics we will specifically consider parameter values obtained from fits at 



one-loop order, which are available in the literature [8j, ll2|, [13j . We begin by noticing that 
our results in Eqs. (jSJ), ffl3]) . and f|T5|) show that / is a common factor affecting the overall 
normalization only. Similarly, C is a common factor, except for the first term in Eq. (I15dl) . 
which is numerically small. Consequently, we fix / and C to their tree-level values, noting 
that the resulting decay rate scales with an overall factor C 2 /f A . In addition, we keep 7 g at 
its value in Eq. (fT2l . as it is well determined. Thus, the ranges of the strong parameters we 
consider are 



0.21 < D-F < 0.34, -2.4 <H< -1.6. (17) 

On the other hand, since the range of the weak parameter he from one-loop fits is large [l3| , 
-2 < 10 7 h c % 4, we let it vary so as to reproduce the experimental decay rates. 

In Fig. E(a) we display the branching ratios calculated from the leading-order (LO) and 
NLO amplitudes above. The black (dark gray) band in the figure shows the effects of 
the parametric variations given in Eq. (|T7|) on the branching ratio obtained from the LO 
amplitude alone (the LO amplitude and only the 7 8 terms in the NLO amplitude). The 
light-gray region results from the LO and NLO amplitudes considered above and varying 
the parameters according to Eqs. ffTBl) and ffTTl) . The dotted lines in this figure bound 
the range 3.3 < 10 4 i3(f2~ — > H~7r + 7r~) < 3.9 implied by the preliminary HyperCP data. 
Evidently, this data can be reproduced in the three cases. 

The corresponding m 3 - 7r + distributions are plotted in Figs. Mb) and (c) for h c < and 
he > 0, respectively, with the variations of the other parameters for the different bands 
being the same as in Fig. Ufa). The he ranges used in (b) and (c) are 0.84 < 10 8 \hc\ < 0.92 
for the black bands, -1.05 < 10 s h c < -0.90 and 0.55 < 10 8 h c < 0.65 for the dark-gray 
bands, and —1.8 < 10 8 he < and < 10 8 he < 1.4 for the light-gray bands, all of which 
have been inferred from the corresponding bands in (a). The figures indicate that some 
softening of the S* dominance in the spectrum is possible with the inclusion of higher-order 
contributions. 
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FIG. 5: (a) Branching ratios for f2 _ — > H~7r + 7r~ and (b,c) the corresponding distributions of 
H~7r + invariant-mass. The black (dark gray) bands come from the LO amplitude only (the LO 
amplitude and the 7 8 terms in the NLO amplitude), and the light-gray bands result from the LO 
and NLO amplitudes we consider, as described in the text. The dotted lines in (a) bound the range 
implied by the preliminary HyperCP data. The dashed curves in (b) and (c) have been reproduced 
from Fig. [3j 

IV. CONCLUSIONS 

We have evaluated the decay Q~ — > H~7r + 7r~ in heavy-baryon chiral perturbation theory. 
At leading order, we found a spectrum dominated by the S*(1530), as had been suggested 
before. This shape is in conflict with the recent preliminary data from HyperCP. The total 
branching ratio is also in conflict with experiment for the central values of C and he, but it 
suffers from a large parametric uncertainty. This uncertainty, however, does not affect the 
shape of the rri'=- 7T + invariant mass distribution. 

A complete calculation at next-to-leading-order contains too many unknown parameters 
to be phenomenologically useful. We have investigated the effect of the NLO corrections 
in three different ways. First, we considered the diagrams in which the weak transition 
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occurs in the meson sector. These corrections are induced by the low-energy constant 7 8 
which is known from kaon decay. Second, we considered the NLO terms in the weak chiral 
Lagrangian which introduce three new effective constants. We studied the effect of these 
constants by varying their value between zero and the value suggested by naive dimensional 
analysis. Third and last, we varied the LO parameters in ranges that included their values as 
determined from tree-level and one-loop fits to other hyperon decay modes. The difference 
between the two kinds of fit is indicative of the size of NLO counterterms that we have not 
included explicitly. When all these factors are considered, we have found that it is possible 
to lower the branching ratio and soften the importance of the H* in the 7713-.^+ distribution, 
as suggested by the data. Beyond this, we can only encourage the HyperCP collaboration 
to fit their data to our result, given in Eqs. (ED, ( 1T31) . and ( 1T51) . 
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APPENDIX A: DERIVATION OF NEXT-TO-LEADING-ORDER WEAK LA- 
GRANGIAN 



The NLO weak Lagrangian generating the weak Vt~Ttp vertices generally contains the 
Dirac structures T M v ■ AT a and T> M 2S ■ AT^. The only possible SU(3) building-blocks 
needed to construct it are therefore the tensors T a b c , Ad e , and Tf g h- Employing standard 



techniques [l4|, we treat the combination T abc Ad e Tf g h as a tensor product (10 <g> 8) ® 10. 
Thus we find four different operators that transform as octets, whose irreducible represen- 
tations are 



where 



/ /r\ \ — rp ( /r\ \ 1 rp 1 c 1 rp 

\ l) ab ^bmn®m,r anr ' \ 2J a0 bmn amn 3 ab mno mno ' 

(O ) = f T T (O ) =8 T 

\ "S)ab bmn ' am, op nop > V a,bmno mno 5 



(Al) 



®a,b ^amn ^bmo ^on ' ^abc ^abm ^mc ^acm ^mb ^bcm ^ma ' (^^) 

^~ab,cd ^cdm i^amo ^bo ^bmo *^ao) 

5 i^aty^dmn ^bmo ^ adJ^cmn ^bmo ^bc^dmn ^amo ^bdr^cmn ^amo) ^no ' (^^) 

@a,bcde = ^bcd ^ae + ^bce ^ad + ^bde ^ac + ^cde ^ab 

— ~fi,(T A 4-T A 4-T A ) 

6 ab V cdm "Trie ' cem "^md ' dem "^mc ) 

— -S(T A 4-T A 4-T A } 

6 u ac \ bdm "^me ' bem "^md ' dem "^mb ) 

— -S At a a-T a 4-T a ) 

6 ad V bcm "me ' cem "^mb ' bem ^-mcj 

6 ^ae {^bcm ~^md ^cdm ^mb ^bdm -^mc) • (A4) 



10 



The tensors (0i,2,3,4)a& and o ab are all traceless, d abc is fully symmetric in its indices, f abfid 
satisfies the symmetry relation f abtCd = Ua,cd = ^ab,dc = fba,dc an d tracelessness condition 
i~ab,cb — 0, and 6 at bcde is symmetric in its bcde indices and satisfies 6 a ,abcd — 0. 

The only possible building blocks needed to construct the NLO weak Lagrangian gener- 
ating the weak Q~Bip vertices are the tensors B ab , A cd , and T de f. Treating the combination 
B a b Acd Tdef as a tensor product (8 <g>8) ® 10, we find four different operators that transform 
as octets, whose irreducible representations are 

(^l)ab ^bmn^mo^ano ' (^ 2 ) a fe ^bmn^'moT ano , 

—J t _ I A f t ((0'\ = f T T 

\ 3/ ab bmn amn ■} ab mno mno ' V 4/a6 bmn am,op n 



(AS) 



■ nop , 



where 



V* = {B,A} ab -l(BA)5 ab , T ab = [B,A] ab , (A6) 

tabc ~ e amn {B mb A nc B mc A nb j + £ bmn [B m( A na + B ma A n( ^ 

+ e cmn {B ma A nh + B mb A na ) , (A7) 
^&,cd = B ac A bd + B ad A bc + -B 6c .A ad + B bd A ac 

bd + &ad&bc) {BA) . (A8) 

The tensors (^i 2,3,4)06, ^a&, and jF ab are all traceless, t a b c is fully symmetric in its indices, 
and f abM satisfy the symmetry relation f abM = T ba>cd = T ab4c = T ba>dc and tracelessness 
condition T ab>cb = 0. 

The resulting NLO weak Lagrangian that contributes to fi~ — > En and transforms 
as (8 L , 1r) is then 

C = (^hah 1 1 + h 2 2 + h 3 3 + h A 4 )) 

+ (\ 0i + h 2 6 2 + fe 3 o 3 + fe 4 o 4 ) } 

+ + K 0' 2 + h' 3 O f 3 + ti 4 0' 4 )) , (A9) 

where Oi and Oi contain the Dirac structures v ■ AT^ and T^2S ■ AT^, respectively, 
and h, L , and h\ and are free parameters. Expanding the Lagrangian yields 

+ ^d^ + E°n- + ••• , (A10) 



where 



(~) (~) (~) (~) 

(~) hj_ _/z_2 _ ' h 3 _ hj 

/»na.»- 5v /g 2v ^' 

3 5 



(All) 
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